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Abstract. From the relic density measurement by WMAP the WIMP annihilation cross section 
can be determined in a model independent way. If the WIMPS are postulated to be the neutralinos 
of Supersymmetry, then only a limited region of the supersymmetric parameter space matches this 
annihilation cross section. It is shown that the resulting positrons, antiprotons and gamma rays 
from the neutrahno annihilation (mainly into bb quark pairs) provide the correct shape and order 
of magnitude for the missing gamma and hard positron fluxes in the Galactic Models and are 
consistent with antiproton production. 
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1 Introduction 

Cold Dark Matter (CDM) makes up 23% of the en- 
ergy of the universe, as deduced from the temperature 
anisotropics in the Cosmic Microwave Background in 
combination with data on the Hubble expansion and 
the density fluctuations in the universe PP . The nature 
of the CDM is unknown, but one of the most popular 
explanation for it is the neutralino, a stable neutral 
particle predicted by Supersymmetry |21|3]. The neu- 
tralinos are spin 1/2 Majorana particles, which can 
annihilate into pairs of Standard Model (SM) parti- 
cles. The stable decay and fragmentation products are 
neutrinos, photons, protons, antiprotons, electrons and 
positrons. From these, the protons and electrons are 
drown in the many matter particles in the universe, 
but the antimatter may be detectable above the back- 
ground from nuclear interactions, especially because of 
the harder positron and gamma spectra expected from 
neutralino annihilation. This so-called indirect detec- 
tion of Dark Matter has been discussed much before 
(see e.g. Ref. j^). Our results differ from these previous 
results by performing a statistical analysis to gamma 
rays, antiprotons and gamma rays simultaneously and 
taking into account the best known propagation mo- 
dels and all constraints from WMAP and electroweak 
data on the SUSY parameter space. More details of 
this analysis can be found in the contributed paper to 
this conference i5|. 



2 Annihilation Cross section Constraints 
from WMAP 

In the early universe all particles were produced abun- 
dantly and were in thermal equilibrium through an- 
nihilation and production processes. At temperatures 



below the mass of the neutralinos the number density 
drops exponentially. The annihilation rate F =< av > 
n-^ drops exponentially as well, and if it drops below 
the expansion rate, the neutralinos cease to annihilate 
and a relic cosmic abundance remains. For the case 
that < crw > is energy independent, which is a good 
approximation in case there is no coannihilation, the 
present mass density in units of the critical density is 
given by 
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One observes that the present relic density is inversely 
proportional to the annihilation cross section at the 
time of freeze out, a result independent of the neu- 
tralino mass (except for logarithmic corrections). For 
the present value of H^h^ = 0.11 the thermally aver- 
aged total cross section at the freeze-out temperature 



of must have been 3-10 



This can be 



achieved only for restricted regions of parameter space 
in the MSSM, as will be discussed in the next section. 
Note that the annihilation cross section is given by the 
Hubble expansion and therefore not dependent on the 
WIMP model. 



3 Predictions from Supersymmetry 

The mSUGRA model, i.e. the Minimal Supersymmet- 
ric Standard Model (MSSM) with supergravity inspired 
breaking terms, is characterized by only 5 parameters: 
iTT-Oi "^1/2: tan/3, sign(/i), AojEI- Here mp and mi/2 
are the common masses for the gauginos and scalars 
at the GUT scale, which is determined by the unifica- 
tion of the gauge couplings. Exact gauge unification is 
still possible with the precisely measured couplings at 
LEP 0. 
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Fig. 1: The light shaded (blue) area is the region allowed 
by WMAP and the contours of larger f2h^ are indicated 
by the dashed lines in steps of 0.05. The upper plot is 
for tan j3=51 and Aq — mo, while the lower plot is for 
tan /3 =55 and Aq — 0. For the last parameters the neu- 
tralino annihilation hits the pseudoscalar Higgs resonance, 
which allows heavier neutralinos with still a small enough 
relic density. The black dots indicate the resonance region, 
where |mA — 2m;^o| < 10 GeV. The excluded regions, where 
the stau would be the LSP or EWSB fails or the boost fac- 
tors are above 10 are indicated by the dots. 



The neutralinos, which are assumed to be the sta- 
ble, lightest supersymmetric particles, can annihilate 
through higgs- and Z-exchange in the s-channel and 
SUSY particles (neutralinos, charginos, sfermions) in 
the t-channel. At large values of tan [3 the dominant 
channel is the pseudoscalar Higgs exchange with bb 
quarks in the final state, which lead to a well defined 
shape of the final state gammas, positrons and elec- 
trons, since the annhihilation is practically at rest. 
The regions of parameter space allowed by the WMAP 
data are plotted in Fig. ^ for two values of tan (3. It 
is clear that for tan/3 w 50 only a small region is al- 
lowed. Scanning over all possible values of tan/3 the 
neutralino masses allowed by the WMAP data and 
electroweak constraints are in the range of 150-400 
GeV|3], if we exclude the coannihilation regions, which 
would lead to anomously large boost factors, as dis- 
cussed in the next section. For the fits discussed below 



we use a typical mass of 200 GeV, which corresponds 
to mi/2 « 500 GeV. The data are not yet sensitive 
enough to distinguish between masses in the range 
given above. 

4 Global Fits to positrons, antiprotons and 
gamma rays 

Trying to disentangle the contributions from nuclear 
interactions and neutralino annihilation to the anti- 
matter fiuxes and gamma rays is in practice not easy. 
We use the following strategy: the shape of the back- 
ground is taken from the GALPROP program, which 
represents a detailed simulation of our galaxy j8j. The 
main background of hard gammas comes from ttq de- 
cays, which are produced in nuclear interactions and 
inverse Compton scattering of electrons on photons. 
The shape of the neutralino annihilation signal is taken 
from DarkSusyl^. These shapes are then multiplied by 
an arbitrary normalization factor, which is left as a free 
parameter in the fit to the data. 

The following data were used in the fit: 

— Gamma ray data from the galactic center in the 
angular range 330° < £ < 30° and -5° < 6 < 
5° from the EGRET space telescope, which has 
been taking data for about 9 years on the NASA 
Compton Gamma Ray Observatory (CGRO). We 
use the data as presented in Ref. 10 . 

— Positron data from AMS pi and HEAT [Hj. 

— Antiproton data from BESS in the years 1997 and 
1998 iJJj 

The fit results are shown in Fig. [21 The free param- 
eters are only the normalization factors for signal and 
background for each of the particle species and their 
values have been indicated in the figures. The boost 
factors, i.e. the free normalization factor after correct- 
ing for the different propagations and energy losses, for 
antiprotons, positrons and gamma rays are all around 
5-7 for the NEW halo profile 14 taken^ Much larger 
factors are not expected from theories of galaxy for- 
mation. If we select SUSY parameters in the so-called 
coannihilation region, where e.g. the stau and neu- 
tralino are almost degenerate, the boost factors come 
out to be much larger, since the fast annihilation cross 
section in the early universe by stau-neutralino coanni- 
hilation does not operate in the present universe any- 
more and the small present annihilation cross section 
for heavy neutralinos needs a large boost factor to fit 
the data. The regions for which the boostfactors are 
above 10 are indicated in Fig. ^ 

The improves significantly with the inclusion 
of Dark Matter in the fits. The y^/d.o.f. is reduced 
from 113/35 (110/38)for the background only fit to 
29/32 (33/35) for the fit including neutralino annihi- 
lation, where the numbers in brackets are valid, if one 
takes the shape and normalization from GALPROP, 
while the first numbers are obtained if only the shape 
is taken and the background normalization is left free. 

^ We use the default (a,/3, 7) =(1,3,1) for a scale a = 10 
kpc and a local relic density of 0.8 GeV/cm^. 
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Fig. 2: From top to bottom: Gamma ray , positron and an- 
tiproton spectrum with contributions from nuclear interac- 
tions (grey/yellow) and neutralino annihilation (dark/red) 
for a neutralino mass of 207 Ge V. The normalization fac- 
tors for signal, called boostfactor, and background (bg scal- 
ing) and the values of x^) with and without signal have 
been indicated. 



This corresponds to about a 4 (6) cr effect, if calculated 
with Gaussian errors. For the antiprotons the increase 
in probability is the least significant, as expected, since 
the shape of background and signal are similar. 

It should be noted that the statistical significance 
is independent of the choice of halo or propagation pa- 



rameters, since different halo or propagation parame- 
ters would only lead to different normalization factors 
in the fit, but the is not affected, since it is only 
sensitive to the shape of the distribution with free nor- 
malization parameters. 



5 Conclusion 

It is shown that the discrepancies between EGRET 
data and the galactic models can be reduced by taking 
as an additional source of hard gammas the annihila- 
tion of Dark Matter, assuming Dark Matter is made 
of neutralinos, as predicted by Supersymmetry. In ad- 
dition, it is shown that adding the positrons from neu- 
tralino annihilation in the same Dark Matter model to 
the same background model improves also the fit 
to the positron data significantly, while the increase in 
antiprotons is compatible with the data. These facts, 
statistical significant improvement of the global fit for 
positrons, antiprotons and gamma rays simultaneously 
for a supersymmetric model with an annihilation cross 
section compatible with the model-independent WMAP 
value, provide strong experimental evidence for the su- 
persymmetric nature of Dark Matter. 

We thank V. Moskalenko and A. Strong for sha- 
ring with us all their knowledge about our galaxy, O. 
Reimers to provide us with the EGRET data. This 
work was supported by the DLR (Deutsches Zentrum 
fiir Luft- und Raumfahrt). 
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